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Abstract. We deal with operators in R" of the form 




where a(jt), b(x) are positive, bounded and periodic functions. We denote by L per the set of such operators. 
The main result of this work is as follows: for an arbitrary L > and for arbitrary pairwise disjoint intervals 
(aj,f3j) c [0, L], j = 1, .. .,m (m e N) we construct the family of operators {A e € L per } e such that the 
spectrum of A £ has exactly m gaps in [0, L] when £ is small enough, and these gaps tend to the intervals 
(aj,(3j) as e — > 0. The idea how to construct the family |A e } £ is based on methods of the homogenization 
theory. 

Keywords: periodic elliptic operators, spectrum, gaps, homogenization. 



Introduction 

Our research is inspired by the following well-known result of Y. Colin de Verdiere flU: for 
arbitrary numbers = A\ < A 2 < ■ ■ ■ < A m (m e N) and n e N \ {1} there is a n-dimensional 
compact Riemannian manifold M such that the first m eigenvalues of the corresponding Laplace- 
Beltrami operator -A M are exactly Ai,..., A m . In the work | fl6| we obtained an analogue of this 
fact for non-compact periodic manifolds: for an arbitrary m pairwise disjoint finite intervals on the 
positive semi-axis (m e N) a periodic Riemannian manifold is constructed such that the spectrum 
of the corresponding Laplace-Beltrami operator has at least m gaps, moreover the first m gaps are 
close (in some natural sense) to these preassigned intervals. 

The goal of the present work is to solve a similar problem for the following operators in W 
(n > 2): 

1 " f> I f> \ 
A = -b- ! div(aV) = -— V — a(x)— , a,b e H per 

b(x) dx k \ dx k j 

where H per is a set of measurable real functions in R" satisfying the conditions 

{3C~, C + > : C < f(x) < C + , Vx e R" (boundedness from above and form below) 
Vi G Z", Vjc e R" : f(x + i) = f(jc) (periodicity) 

The operator A acts in the space L 2; b(R n ) = lu e L 2 (R"), ll"lli 2b(Rn) = / \u(x)\ 2 b(x)dx >, it is self- 

adjoint and positive. We denote by L per the set of such operators. 

l 



f € H per 



Operators of this type occur in various areas of physics, for example in the case n = 3 the 
operator A governs the propagation of acoustic waves in a medium with periodically varying mass 
density (a(x)) -1 and compressibility b(;c). 

It is well-known (see e.g. [[T71 ) that the spectrum <x(A) of the operator A e L per has band struc- 
ture, i.e. cr(A) is the union of compact intervals [a k ,a k ] c [0, oo) called bands (a~ = 0,a k /" oo). 

k— >oo 

In general the bands may overlap. The open interval (a,fi) is called a gap if (a,/3) n <x(A) = and 
a,p E cr(A). 

The main result of this work is the following 

Theorem 0.1 (Main Theorem). Let L> Obean arbitrary number and let (aj,fij) (j= 1 , . . . , m, m e 
NJ oe arbitrary intervals satisfying 



< »i, Qfj < < j = 1, m - 1, or m < fi m < L (0.1) 

Lef n e N\ {1}. 

77zen one can construct the family of functions ja £ e H per } an J the function b 6 H per swc/z f/ia? 
the spectrum of the operator A £ = b _1 div(a £ V) has the following structure in the interval [0, L] 
when s is small enough: 

t m \ 



o-(A £ )n[0,L] = [0,L]\ 



(0.2) 



where the intervals (a ,£ ,/?p satisfy 



Vj = l,...,m: lima £ = a y , lim/? e =/?,• (0.3) 

e— *0 ^ £— >0 J 

Moreover, a E (x), b(x) are step-functions having at most m + 1 values. 

Remark 0.1. It follows from (10. 1 b - (IQ.3b that the operator A £ has exactly m gaps in [0, L] when 
s is small enough. In general, the existence of gaps in the spectra of operators from L per is not 
guaranteed, for instance in the case of constant a(x), b(x) the spectrum o~(A) coincides with [0, oo). 
Various operators from L per with gaps in their spectrum were studied in the works ll5]- n"Tll2~2ll27Tl 
(see also the overview IfTllO . In these works spectral gaps are the result of high contrast either in 
the coefficient a(*) E|91[ni|22 or in the coefficient b(*) EE) or in both coefficients ll5l[T0ll22l 
(the last three works deal with the Laplace-Beltrami operator in W with conformally flat periodic 
metric; obviously, this operator belongs to L per ). 

The operator A £ constructed in the present work also has high contrast in the coefficients (namely, 

lim ( — 3 — | = oo), but their form essentially differs from the form of the coefficients in the 

e^o \ mm xeR „ a £ (x) / 

works mentioned above. 

The idea how to construct the functions a £ (x), b(x) has come from the homogenization theory. 
We briefly describe this construction. 

m 

Let s > be a small number. Let G £ = |J |J G £ . be a union of pairwise disjoint spherical shells 

i'eZ" j=\ 

Gf. lying in R". It is supposed that the following conditions hold (see also Fig. 1): 

• for any fixed j e {1, . . . ,m} the shells G £ . are centered at the nodes of e-periodic lattice in 

R", 

• the shells G e Qj (j = 1, . . . , m) belong to the cube {x = (x\, . . . , x n ) el": < x k < e, Vfc}. 



Fig. 1. 



The external radius of the shells is equal to r E = rs (r > 0), the thickness of their walls is equal to 

m 

d E = s y (y > 3). By B E . we denote the sphere interior to G?.. We set B £ = \J \J B E .. 

i€Z" 7=1 

We define the functions a e {x), b e {x) by the formulae 

a E (x) = l 1, xt^"\G E , b s (x) = [l> xgR"\(B £ UG e ), 

\a^\ xeGf, W [bj, xeBfjUGfj, 

where aj, bj (j = 1, . . . , m) are positive constants, which will be chosen later on. We consider the 
operator 

3? = -(^T'div^V) = VfUd] 

It will be proved (see Theorem 11.1 I below) that the spectrum of 3K £ converges to the spectrum of 
some operator 31° acting in the Hilbert space L 2 (R") ©_ Z^/o-^R"), where pj, crj (j = 1, . . . ,m) 

j=l,m 

are positive constants. The spectrum of 31° coincides with the set [0, 00) \ ( \J (crj, fiM, where the 
intervals (crj, fij) satisfy 



0<cr u crj < jj.j < o-j+u j = l,m- 1, cr m </J, m < 



00 



and depend in a special way on cij and bj. 

More precisely, we will prove that for an arbitrary L > fi k the spectrum of the operator 3K £ has 
the following structure in the interval [0, L] when s is small enough: 

( m \ 



<r(3l E )n[0,L] = [0,L]\ 



where the intervals (o^,^) satisfy 

Vj = 1 , . . . , m : lim cr E = cr,-, lim /i £ = /i ,■ 

e— »0 e— >0 J 



Furthermore, we will prove (see Theorem 1 1 .21 below) that for arbitrary intervals (aj,/3j) (j = 



1, . . . ,m, m e N) satisfying (10.11) one can choose such ai, bj in (10.41) that the following equalities 
hold: 

Vj = 1, . . . , m : o-j = aj, Hj = fij (0.5) 

Finally we set (below y e R' ! ) 

a £ (y) = s~ 2 (f{x), b(y) = b E {x), where x = ve 

(obviously, b(y) is independent of s). It is clear that a e , b belong to H per and are step-functions 
having at most m + 1 values. It is easy to see that the spectra of the operator 

A e = b _1 div(a e V) 

and the operator J?l e coincide (in fact, A £ is obtained from Jl E via change of variables x = ys). 
It follows from Theorem OP that cr(A E ) satisfies (lO2l)-(l03T). 

We remark that the gaps open up in the spectrum of A e because of the high contrast in the 
coefficient a e (x). The coefficient b(x) is independent of s and it is needed only in order to control 
the behavior of the gaps as s — » 0. In fact, the operator -div(a £ V) also has at least m gaps when s 
is small enough, but in general they do not converge to (<*/,/?/) as e — » 0. 

Heuristic arguments. The classical problem of the homogenization theory (see e.g. Ifn-T51 [T%1I2"4T - 
l26ll ) is to describe the asymptotic behaviour as e — » of the operator Jl E which acts in L 2 (Q) 
(fl c R" is a bounded domain) and is defined by the operation 

^ = -div (a e V) 

and either Dirichlet or Neumann boundary conditions on dQ. Here 

a s (x) = a(xe _1 ), where a e H per (0.6) 

It is well-known that Jl E strongly resolvent converges to the operator (so-called "homogenized 
operator") 



dxtdxj 



k,i=\ 

where the constants a*' satisfy: 3C , C + > s.t. V£ e R" C^l^l 2 <fl H f Jt f, < C + |^| 2 . 

It is interesting to study the asymptotic behaviour of the operator Jl E when a £ has more com- 
plicated form comparing with (|0.6I) . In particular interest is the case when a E is bounded below 
but not uniformly in s. This is just our situation (see (10.41) ): for fixed e one has mma E (x) > 0, 

but lim(mina e (x)] = 0. Such type problems were widely studied in lfT8l Chapter 7]. In particu- 



+0 \j 



lar, the authors considered the operator jR^ E which acts in L 2 (Q) and is defined by the operation 
J{®' £ = -div(a e V) and the Dirichlet boundary conditions on d£l. Here Q, c R" is a bounded 
domain, a E is defined by (10.41) (only the case m = 1 was considered). It was proved that J?t^' £ 
converges as e — » (in some sense which is close to strong resolvent convergence) to the operator 
acting in the space L 2 (Q) © L 2 , p /o-(Q.) and being defined by the operation 

<°=("^; p i\ (0.7) 

and the definitional domain D{^f) = {(w, v) 6 H 2 (Q.) © L 2i p/o-(Q) : u\aa = 0|. Here ~a,p, cr are 
positive constants that do not depend on Q. A similar result is valid for the operator Jl^ £ (the 



5 



superscripts "D" and "N" mean Dirichlet and Neumann boundary conditions): the corresponding 
homogenized operator J?T^'° is defined by operation (10.71 ) and the definitional domain D(J?l^' ) = 

{(u, v) e H 2 (Q) e L 2 , pAr (£l) : || an = 0). 

Although in general the strong resolvent convergence of operators does not imply the Hausdorff 
convergence of their spectra (see the definition at the beginning of Section [5]), but suppose for a 
moment that this is true for the operators Jl^ £ and tft^f , i.eQ 

cr(J{°' £ ) cr(J{° fl ), cr(Ji^ E ) -> cr(^ ) in the Hausdorff sense 

We denote fi s = ji € R" : \x\ < R). One can prove (for example, it follows from lfl"5l Proposition 
2.3]) that 

V£lcR": (cr, yt/) n cr(Jl° /Nfi ) = 

V[<T, J + ] c [0, oo) \ (cr, jS) 3R d >0 : a(3& Nfi ) n [<T, <T] * for i? > 

where D/N is either D or N, \x = cr + p. These suggest that when e is small enough the operator 
Jl E has a gap in the spectrum and this gap tends to the interval (cr,//) as e — > 0. 

The close problem was also considered in ETI where the authors studied the asymptotic be- 
haviour of the attractors for semilinear hyperbolic equation d 2 tt u + tfl^f'u + f £ (u) = h £ . 

We remark that the proof of the resolvent convergence in lfT8ll is based on the method of so-called 
"local energy characteristics". This method is well adapted for both periodic and non-periodic 
operators but it is quite cumbersome. Therefore in the present work following lfT6l we carry out 
the proof in more simple fashion via the substitution of a suitable test function into the variational 
formulation of the spectral problem. 

In the next section we describe precisely the operator Jl E and formulate Theorems ll.im.2[ Their 
proofs are carried out in Sections [2HTJ 

1. Construction of operators Jl E and main results 

Let neM\{l),meN. Let the points Xj e R" (j = 1, . . . , m) and the number r > be such that 
the closed balls Bj = jx s W : \x - xj\ < rj are pairwise disjoint and belong to the open cube 

Y = {x = (x u . . . , x n ) £ R." : < x k < 1, Vfc} 
Let s > 0. We introduce the following notations (below i e Z", j = 1, . . . , m): 

x £ j = s(xj + i) 

G E j = [x 6 R" : r £ -d E <\x- x E j\ < r E ) , B E j = [x 6 R" : \x - *f y | < r E - d E ) 

where 

r E = rs, d E = s y , y > 3 

We also denote 

m m 

G£ = UU G ^ fiC = UU 5 'i' f e = ^\{g^JW) 

ieZ" y=l ieZ" ;=1 



1 We will prove this statement in Section [5] (the only difference is that we will consider quasi-periodic boundary 
conditions, but for Dirichlet and Neumann boundary conditions the proof is similar.) 
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We define the piecewise constant functions a £ (x), b £ (x) by the formulae 

fl, xeF £ UB £ , 
a £ (x) = \ ^ (1.1) 

[a £ = ajS y+ \ x 6 G £ j (i 6 IT, j = 1, . . . ,m), 

fl, xeF £ , 

b £ (x) = \ (1.2) 
[bj, xeBfjUG*. (ieZ", j= l,...,m), 

where a ; -, bj (j = 1, . . . , m) are positive constants. 

Now we define precisely the operator By L 2 ^(K") we denote the Hilbert space of functions 
from L 2 (W l ) with the following scalar product: 

(u, v)l 2>a£: (r») = J u(x)v(x)b £ (x)dx, 

Remark that 

C || • ||l 2 (R") ^ II ' IIl 2 ,, £ (R") < C + || • ||l 2 (r«) (1-3) 

where the positive constants C , C + are independent of s. By Tjt n [u, v] we denote the sesquilinear 
form in L 2 foE (R") which is defined by the formula 



Ve." v] = ^ a £ (x) (Vw, Vv) dx 



, _ v-i du dv 
with dom{ri £ ) = H (W). Here (Vw, Vv) = > . The form is densely defined, closed 

dx k dx k 

and positive. Then (see e.g. Q3)) there exists the unique self-adjoint and positive operator Jl £ 
associated with the form tj^„ [u, v], i.e. 

(W £ u, v) L2 ^( R n) = r] £ ,,[u, v], Vw e dom(j?l e ), Vv e dom^,,) (1.4) 

Its domain dom(J?l e ) consists of functions u belonging to the spaces H 2 (F £ ), H 2 (G £ j), H 2 (B £ p (for 
any i e Z", j = 1, . . . ,m) and satisfying the following conditions on the boundaries of the shells 

(ur = (uT and *^(bJuG*), ^ 

(ur = (u)- and = (|) , xe^ 



where by + (resp. -) we denote the values of the function u and its normal derivative on the 
exterior {resp. interior) side of either d \B S .. U Gf,) or dB £ ... For sufficiently smooth w the operator 
J?l £ is defined locally by the formula 

By cr(Jl £ ) we denote the spectrum of the operator M E . In order to describe the behaviour of 
cr(Jl £ ) as s — > we introduce some additional notations. 



In the domain F = Y \ \J Bj we consider the following problem (below k = l,...,n): 

7=1 



Av^ = 0, x e F 
dv k 

— = n k , x e d 
on 



v k ,Dv k are y-periodic, i.e. Va = \,n : 



for x — (xj , ^2, . . . , 0, . . . , Xn) 



d-th place 



(1.7) 



where n = (m, . . . ,n„) is the outward unit normal to IJ B;, e a = (0, 0, . . . , 1, . . . , 0). It is known 



;'=i 



ot-th place 



(see e.g. (3J) that the unique (up to a constant) solution v k {x) of this problem exists. We denote 
a** = -i- J" (V(x* - v fe ), V(x z - v/)) <ix, fc, / = 1, . . . , n 



The matrix A = j*r J is symmetric and positively defined (see e.g. (31 Chapter 1, Proposition 2.6]). 

Remark 1.1. In the case when m = 1 and the center of ball Bi coincides with the center of the cube 
Y the matrix A = ja*'} has more simple form, namely A = al where I is the identity matrix, H > 0. 
This follows easily from the symmetry of the domain F. 



We denote 



na; 



cr, 



Pj 



aj \dBj\ 



(1-8) 



; rb/ rj \F\ 

We assume that the numbers aj and bj in (ll.ll) - (ll.2l) are such that cr t + crj if i ^ j. For definiteness 
we suppose that crj < o" i+1 , j = \,...,n - 1. 

And finally let us consider the following equation (with unknown A e C): 



T(A) = 1 + V = 

^— l CT — A 



7=1 



(1.9) 



It is easy to prove (see Section H]) that this equation has exactly m roots fij (j = 1, . . . , m), they are 
real, moreover they interlace with o~j, i.e. 



(Tj < fij < o-j+u j = 1, m - 1, cr„, < ix m < oo 
Now we are able to formulate the theorem describing the behaviour of cr(Jl E ) as s — » 0. 

Theorem 1.1. Le/; L Z?e an arbitrary number such that L > fx m . Then the spectrum o~(j?l £ ) of the 
operator Ji £ has the following structure in [0, L] when s is small enough: 



o-m £ )n[0,L] = [0,L]\ 



W=i 



where the intervals (o^,^) satisfy 



V/ = 1 , . . . , m : lim cr E = cr,, lim = [i ,• 

£— >0 £— >0 ^ 



(1.10) 



(1.11) 



The set [0, oo) \ I Q (cry, ju 7 ) I coincides with the spectrum <x(J?T ) of the self-adjoint operator J{° 
which acts in the space LaiW 1 ) ©_ Z^/o-^R") and is defined by the formula 

j=\,m 



JLa dx k dx, j=1 

cr\(u\ - u) 
cr 2 («2 - u) 

0~m(Mm ~ M ) 



, u 



( u \ 

U\ 

u 2 



e dom(J{ ) = H 2 (R n ) (^L 2 , Pj/(Tj (R n ) 

j=\,m 



To complete the proof of Theorem lO.il we have to choose such aj and bj in (11.11) . (11.21) that (10.51) 
holds. 

Theorem 1.2. Let (aj,Bj) (j = I, ... ,m) be arbitrary intervals satisfying ( 10.71) . 
Then ( 10.51) /io/ds ifwe choose 

\F\ 



aj = 



\oBj\ }_ L \aj-ajj r\dBj\ aj LJ \aj-aj) 



i=\,m\i*j 



Remark 1.2. Since the intervals (aj,8j) satisfy (10.11) then 

V/ : Bj > aj, Vi + j : signOS,- - aj) = signfe - aj) ± 
Bi-af 



Therefore {Bj - aj) f] 

i=l,m\fc j 



at - a. 



> and thus the choice of aj and bj is correct. 



The scheme of the proof of these theorems is as follows. 

In Section[2]we introduce the functional spaces and operators that are used throughout the proof. 
Also we present well-known results describing the spectrum of the operator Jl £ 
In Section [3] we prove several technical lemmas. 
In Section H] we show that 



o-m ) = to, oo) \ 



(1.13) 



Section \5\ is a crucial part of the proof: we show that as s — > the set a{^R £ ) converges in the 
Hausdorff sense to the set cr(Jl Q ). 

In Section [6] we prove that for an arbitrary L > the spectrum cr{J{ E ) has at most m gaps within 
the interval [0, L] when s is small enough. Together with the Hausdorff convergence this fact 
implies the statements of Theorem ll.il 

And finally in Section [7] we prove Theorem 1 1.21 

Remark 1.3. We present the proof of Theorem II .11 for the case n > 3 only. For the case n = 2 
the proof is repeated word-by-word with some small modifications (for example in formula (13.101 ) 
below r 2 " has to be replaced by In r). 

2. Preliminaries: functional spaces and operators 

Below Q is a domain in R" with Lipschitz boundary (if dQ. ± 0), for simplicity we suppose that 
d£l n |J G E = 0. Throughout the paper we will use the following functional spaces: 

ij 
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• L 2j ^(Q) be the Hilbert space of functions from L 2 (£l) with the scalar product 

(w, v)l 2 b e(ci) — I u(x)v(x)b E (x)dx 

LI 

o 

• H l (Q.) be the subspace of H l (Q.) consisting of functions vanishing on dd, 

o 

• C°°(£X) be the space of functions from C°°(£X) compactly supported in CI, 



iI 2 ' s (Q.) be the space of functions belonging to H 2 (£l n G?.), H 2 (a. i i yi ^ ^ , j - 
1, . . . , m), // 2 (Q n F E ) and satisfying conditions (11.51) for all shells Gf. belonging to f2, 
C 2£ (Q) be the space of functions belonging to C 2 (Q n Gf.), C 2 (Q n flf.) (i e Z", / = 



1, . . . , m), C 2 (Q Pi i 7£ ) and satisfying conditions (11.51) for all shells G £ . belonging to Q. 
For w, v e H l (Cl) we denote 

?7q[m,v] = J" a £ (x) (Vm, Vv) dx (2.1) 

By ?7^' £ (resp. v^) we denote the sesquilenear form defined by formula (12.11) and the definitional 

domain ^(XX) {resp. H l (Q)). 

Similarly to the operator 3l E (see (11.41) ) we define the operator 3^ E (resp. 3^ £ ) as the operator 
acting in L 2 and associated with the form ?7^' £ (rasp. ^ e ). The definitional domain dom(J?T^' £ ) 
(resp. dom(J?I^' £ )) consists of functions from H 2e (Q.) satisfying the condition §^| 5n = (res/?, 
"bn = 0) that justifies the upper index "N" (resp. "D") which indicates the Neumann (resp. 
Dirichlet) boundary conditions. 

The spectra of the operators Jl^ E , 3^ E are purely discrete. We denote by |/l^' £ (f2)}^ gN (resp. 

{^f - s (Q)} feN ) the sequence of eigenvalues of J?t^' £ (resp. written in the increasing order and 

repeated according to their multiplicity. 

Now let us describe the structure of the spectrum cr(Jl E ) of the operator jft E . The operator Jl E is 
periodic with respect to the periodic cell 

Y E = {x 6 R" : < x k < s, V/fc} 

We denote T" = {9 = (0 lf . . . , 6 n ) e C : \0 k \ = 1, V£}. For e T" we introduce the functional 
space H\(Yq) consisting of functions from H l (Y E } ) that satisfy the following condition on dY^\ 

Vfc = 1, n : w(x + ee^) = 9uu(x) for x = (x\, x 2 , . . . , 0, . . . , x n ) (2.2) 

T 

&-th place 

where e k = (0, 0, . . . , 1 , . . . , 0). 

By r]y E we denote the sesquilenear form defined by formula (12.11) (with Y% instead of Q) and the 

definitional domain H l g (Y E ). 

We define the operator jf' E as the operator acting in L 2 ^(K?) an d associated with the form rj e v E . 
Its definitional domain dom(J?lyf ) consists of the functions from H 2 - E (Yf,) satisfying the condition 





(12.21) and the condition 



du du 

Vfc = 1, n : — (x + se k ) = 6 k —(x) for x = (xi, x 2 , . . . , 0, . . . , x n ) 

&-th place 
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The operator has purely discrete spectrum. We denote by {^(^o)}^ tne sequence of 
eigenvalues of written in the increasing order and repeated according to their multiplicity. 

o 

From the min-max principle (see e.g. ||23~T0 and the enclosure H X (Y^) d H\(Yq) d H 1 (Yq) one 
can easily obtain the inequality 

Vfc 6 N : A* e (Y s ) < A k E (Y E ) < A°' e (Y e Q ) (2.3) 

The following fundamental result (see e.g. IfTTl ) establishes the relationship between the spectra 
of the operators 3l E and jf v f. 

Theorem. One has 

CO 

<r(& s ) = {Jjk(tf £ ) (2.4) 



k=i 



where $ k (Jff~) = U \A\ E {J^\. The sets $ k (J{ E ~) are compact intervals. 

0eT" 

Remark 2.1. It is clear that if e -1 e N then Ji E is also y-periodic operator, i.e. a E (x + i) = a E (x), 
b E (x + i) = b E (x) for any i e IP, x e W. So in this case we have an analogous representation 

oo 

o-m = \JjkW £ ) (2.5) 

k=l 

where J k (Jl E ) = \J {A e k ' E (Y)), A e ^ E (Y) is the k-th eigenvalue of the operator jf Y E which acts in 
L 2 ^(Y) and is defined by the operation (11.61 ) and the definitional domain 

(u(x + e k ) = 9 k u(x) 

du du foTX = (x 1 ,x 2 ,...,0,...,x n ) 

—(x + e k ) = 6 k —(x) t 
OX k OX k fc-th place 

Studying the Hausdorff convergence of cr(Jl E ) as e — > we will use the representation (12.51) . 



while estimating the number of gaps in the interval [0, L] we will use the representation (12.41) . 

3 . Auxiliary lemmas 

In this section we prove some technical lemmas. In order to formulate them we introduce some 
additional notations. 
We denote 



k = — mindist 

2 j=l,m 



( 



Bj, dY U 



\k*j It 

Recall that the closed balls Bj are pairwise disjoint and belong to the open cube Y, hence k > 0. 
We introduce the following sets (below i eZ", j = 1, . . . , m): 

• Yf = {x = (xi, . . . , x n ) 6 R" : is < x k < (i + l)e, Vfc} 

. Ff = yf\u(5pJGj) 

• R E j = jx e I" : r E <\x- xf.| < r 6 + /eg} 



Df. = fx e R" : |jc - xf ; | < r £ + /eel = Bf. U G e . U R 
• S £ . = (jc £ R" : |jc - X s .! = r £ + *e) = <9Z> 



e 
U 



11 



• Cf. = [x g W : |jc - 4I = r' 

• Cf. = [x e R f 



We also denote 
and set 



I E = [i = (h,..., i n ) e Z" : < i k < (e- 1 - 1), V/c) 

m m 
!'eJ e ;'=1 !'eJ e y'=l ;eJ E 

Remark that if e N then Y= \JYf- 

By (u) B we denote the average value of the function u over the domain B c K" (if \B\ + 0), i.e. 



(u) B = — J u{x)dx. If Z c R" is a (n - l)-dimensional surface then the Euclidean metrics in R" 

B 

induces on S the Riemannian metrics and measure. We denote by ds the density of this measure. 

1 r 

Again by {u)z we denote the average value of the function u over S, i.e = — J m^s (here 
PI = / ds). 

Y. 

If v] is a sesquilinear form then we preserve the same notation 77 for the corresponding 
quadratic form, i.e n[u] = n[u, u]. 

By Xa we denote an indicator function of the domain Q, i.e. x a i x ) = 1 for jc e Q and^- n (x) = 
otherwise. 

In what follows by C, Ci... we denote generic constants that do not depend on e. 

Lemma 3.1. Let D be a convex domain in R", d be the diameter of D, X and Y be arbitrary 
measurable subsets ofD. Then for any v e H l (D) the following inequality holds: 

\(v)x-(v)Y\ 2 <C\\Vv\t (D) ^^^ 

Proof. The lemma is proved in a similar way as Lemma 4.9 from lfT8l p.l 17]. □ 

Lemma 3.2. Let s = s N = ± N = 1,2,3 . . . Let v E e H\Y), \\v £ \\ 2 m(Y , < C, v E -» v e H 1 ^) 
strongly in L 2 (Y). Then Vj = l,m: 

y(v £ ) 5 f^ -» -v tfrongfy in L 2 (Y) (3.1) 

y<v E W y£ -> v tfrongfy m L 2 (7) (3.2) 
Proof. For an arbitrary z e J £ and j e {1, . . . , m) one has the following inequalities: 

v£ -< v >llL ( y^ Ce2||Vve|l W) (33) 



s n \(v £ ) Yf - (v £ ) Ff \ 2 < Cs 2 \\Vv% (Yf) (3.4) 

• ij 

s" (v E ) S c. - <v £ )«, <Cs z \\Vv% 2(Rn (3.6) 



< C £ 2 ||Vv £ || / : 2(y£) (3.5) 



2 
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Inequality (13.31) is the Poincare inequality, inequalities (I3.4I) - (I3.5I) follow directly from Lemma [3TTT 
Let us prove inequality (13.61) . We introduce in the spherical coordinates (r, 0), where r is a 
distance to x?., are the angle coordinates. Below by S„_i we denote the (n - l)-dimensional unit 
sphere, by d® we denote the Riemannian measure on S„_i. One has 



v E (r E + ks, 0) - v E (r., 



r°+K£ 

r dv E 

• &>= j -p 



(p, ®)dp, r € (r E , r E + ks) 



We multiply this equality by r" drd®, integrate from r e to r e + ks (with respect to r) and over S„_i 
(with respect to 0), divide by and square. Using the Cauchy inequality we obtain 



kit 



S„_, r"- \ r 



r"- l drd® 



< 



< C 



n 



dv E 

-T-(P,0) 

dp 



p^dpd® 



\ ( r s +K£ 

dp 



r df 
J p"- 



<C^v% im) s 2 - n 



and thus (13.61) is proved. 

It is clear that CD} follows from (Q, (f33]), (EH), and <H2]) follows from (Q, (T3H) . 

Lemma 3.3. The following inequality is valid for an arbitrary v e H l (D E j): 



□ 



lL (GE) <C^- 1 k4v]+^[v] + 



'£2(S*.) 



(3.7) 



Proof. As in the proof of Lemma [3T21 we introduce in G?. the spherical coordinates (r, 0). One has 

r 

—ip,®)dp, re(r E -d E ,r E ) 

Taking into account (11.11) we obtain from (13.81 ) 



(3.8) 



J J \v{r,®)\ 2 r n - l drd® < 2 

S„-i rt-dz 



, J~ r"- l dr ■ (ff- n J |v(r e ,0)| 2 (r £ )"- 1 J0+ 



S n -l 



S„_i V-d" 



dp 



p n ~ l dp 



r dp_ 

J p"- 1 



d® 



< C s r \\v 



.riu,i|2 



'i2(Cf.) 



+ ^- 1 ^,[v]) 



Similarly we obtain 



U)- c (^ 1||v|l W +e||Vv|l W) 



The statement of the lemma follows directly from the last two inequalities. 
Lemma 3.4. lim /if ' e (Df .) = cr,-, where cr, (j = l,...,m) are defined by rii.<SD . 

£->0 J 



□ 
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Proof. Let vf ; e dom(J?I^f) be the eigenfunction corresponding to A"' £ (D £ .) and such that 

'j 



D,£, 



J 

Jb s 



vUx)dx = \B\ 



(3.9) 



Instead of calculating vf in the exact form we construct a convenient approximation for it. 

We introduce in Df. the spherical coordinates (r, 0), r e [0, r £ + /eg). Let ^ : R — > R be a 
twice-continuously differentiable function such that tp(p) = 1 as p < 1 /2 and </?(p) = as p > 1 . 

We define the function v £ . by the formula (below we assume that < r £ - d E that is true for e 
small enough) 



\ re [0,0 

l + A> 2 -"(l-^^g^)), re[ff--d-) 

A £ jr 2 -"<P (^J^) , r e [r £ , r E + ks) 

We choose the coefficients A £ , A £ , A £ , B E in such a way that v £ . satisfies conditions (11.51) : 



(3.10) 



7 1 - a £ L J 



-l 



jji-lgW-l-y 



A £ =A £ 



a £ A £ , 5 £ =A £ (r £ ) 2 -"(a £ -l) 



It is clear that v £ . 6 dom(^f) and J?l £ v £ . = in D £ \ [x : |jc - * £ .| e ^f] U [r £ + f , + ks]}. 
Direct calculations lead to the following asymptotics as s — » 0: 



ll^lkCDM = 0(s n ), 



| v £ - 1 1I 2 + |lv £ II 2 



Using the min-max principle we get 

'/o ^ K/l aj\dBj\ naj 



,£ l|2 



< 



IKyll^^o W rb J 



(Ti 



One has the following estimates for the eigenfunction v £ .: 



s i|2 



(3.11) 
(3.12) 

(3.13) 

(3.14) 
(3.15) 

(3.16) 



The first one is the Friedrichs inequality, the second one is the Poincare inequality and the third 
one follows from Lemma 1331 Furthermore one has the equality 

tfytfj] = ^(^(lIvpL^ + bMtfwp + ^ (llvf. - lH^ + (3.17) 
It follows from (l3TT3l - (|3T7T) that 



rfM] = 0(s n ), 



= o{s n ) as s -» 



(3.18) 



14 



Moreover (l3T9"b . (13.181) imply 



Now let us estimate the difference wf^ = vf. - v £ .. One has 



and thus in view of (13.121) . (13.181) we conclude that 



: .UR E .) 



//'II/. -if/' A' I 



+ llvf { |lL«..^l + 2|||vf,-l|ll,„^ + l|l-v 



£ |,2 1 



\W 



,e l|2 



Furthermore using inequality (13.131) we get 



i£.[wf y ] < -2(^ £ v £ ,vv £ ) L2iE(Df;) + 



v 



,e ll2 



|v e||2 ""W'l^siDfj) 
* y ijHLz^Dfj) 



and in view of (13TTT) . (I3l2l) . (l3TT8T) - (l3T20l) we conclude that 

ij ■> 

The statement of the lemma follows directly from (13.111) . (13.201) . (13.211) . 
Lemma 3.5. lim ^(Df .) = co 

iVoo/ We denote: 

B £ = jy e R" : < \y\ < r - s 7 ~ 1 } , B = {y e R" : < \y\ < r) 

G £ = [y e W : r- s 7 ' 1 < \y\ < r] 

R = {y e W : r < \y\ < r + k) , D = {y e W : < |y| < r + k) 

Also we introduce the functions a £ (y), b(y): 

a £ (y) = a £ (ye + * £ .), b(y) = & £ (ye + x £ .), y e D 

(it is clear that b in independent of s). 

By Ap' £ we denote the operator acting in L 2 ,b(D) and being defined by the operation 



dyj 



and the definitional domain dom(Ap £ ) which consists of functions v belonging to H 2 (B £ ), H 2 (G £ ), 
H 2 (TL) and satisfying the conditions 



(3.19) 



(3.20) 



(3.21) 
□ 



dv 



(vr = (vr and (-] = a £ (-j , y e dB 



dv 



dv 



(vr = (v)- and a £ (-j = (-) , y eeTC 



<9n 
dv 



dn 



v = 0, 



ye 3D 



We denote by /t^' £ (D) the fc-th eigenvalue of the operator Ap £ . It is clear that 

VkN: Af ' £ (D) = e 2 if' £ (D £ .) 



(3.22) 
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Below we will prove that 

Vk 6 N : ^f' e (D) -> A k (3.23) 

where A k is the fc-th eigenvalue of the operator A which acts in the space L 2 (R) © L 2i ^(B) and is 
defined by the formula 



A = - 



A% N 

o b??4 



Here the operator A R ' (resp. Ag) is defined by the operation A and the definitional domain con- 
sisting of functions v £ H 2 (R) (resp. v e H 2 (B)) satisfying the conditions 



V\dD = 0, — 

on 



= (resp. — 
dR\8T> on 



= 0) 

SB 

1 



It is clear that A\ = (Ai coincides with the first eigenvalue of -b ■ : Ag) while 

A 2 > (3.24) 

(/l 2 coincides either with the first eigenvalue of -Ag N or with the second eigenvalue of -Zr^Ag). 
Then the statement of the lemma follows directly from (I3.22I) - (I3.24I) . 

To complete the proof of lemma we have to prove (13.231) . For that we use the following 

Theorem (see [SI). Let 'H £ ,'H be separable Hilbert spaces, let £ £ : <H E -> <H E , £° : <H° -> *K° 
be linear continuous operators, imX° cTc < H°, where <V is a subspace in e H°. 
Suppose that the following conditions C\ - C4 hold: 

C\. The linear bounded operators R £ : Tl — » *H e exist such that ||i? e /ll^/ E — * £ > H/llLo/ or aw .y 
/ 6 < V. Here g > is a constant. 

C 2 . Operators £ £ , £° are positive, compact and self-adjoint. The norms ||X s IU(« e ) are bounded 
uniformly in s. 

C 3 . For any f e <V: \\£ £ R £ f - R £ £°f\\w. -» 0. 

£— >0 

C4. For a/ry sequence f £ e *H £ swc/z ?/zaf sup ||/ e ||^ < 00 the subsequence s' c s and w e *V 
emf smc/? ?/za? ||-£ e / e _ ^wII-h^ — > 0. 

s=e' — >0 



TTzen /or any e 



£— »0 



where {fJ-f}^ and {pk}kLi are tne eigenvalues of the operators £ £ and £°, which are renumbered in 
the increasing order and with account of their multiplicity. 

Let us apply this theorem. We set <H £ = L 2jb (D), «H° = L 2 (R) © L 2 ^(B), £ £ = (A% £ + I)" 1 , 
£° = (A + 1) -1 , <V = K . We introduce the operator R £ : ( H Q -> *H £ by the formula 

[*/](?) = f = (Jk.A)e*° 

[f B (y), yeB, 

Evidently conditions C\ (with q = 1) and C 2 hold. Let us verify condition C3. 
At first we introduce the operator Q £ : H l (B £ ) — » H l (R n ) by the formula 

[<2 e v]O0 = [(2v £ ](^) 

where fc e = (r - e r_1 ) _1 r, the function v £ e /^(B) is defined by the formula v £ (y) = v(y/k £ ) and 
Q : H l (YS) — » //^(R") is the operator with the following properties: 

Vv e ^(B) : [Gv](y) = v(y) for y e B, ||Qv|l ffl(Rn) < CIMI*^ 
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(such an operator exists, see e.g. |[T9l ). One has 

Vv e //' (B £ ) : [Q £ v](y) = v(y) for y e B £ 
Since £ £ ~ 1 as s — » 0, then, obviously, 

Vvg// 1 ^): H2 £ v|| W) < dlMI^ 
Let / = (/ R , / b ) g < K°. We set / £ = R E f, v £ = £ £ f £ . It is clear that 

l|v £ ||L 2 , b (D) < ll/ £ HL 2 , b (D) = ll/ll^o 

One has the following integral equality: 



a £ (y)(Vv £ , Vw £ ) + b(y) (v £ w £ - f £ w £ ) 



dy = 0, Vw £ G H\D) 



Substituting into (13.271) w £ = v £ and taking into account (13.261) we obtain 

a £ \Vv £ \ 2 dy < C 



(3.25) 
(3.26) 

(3.27) 
(3.28) 



Let G // J (R) (resp. Vg G H l (B)) be the restrictions of v £ onto R (resp. the restrictions of 
2 £ v £ onto B). Since v £ G dom(A£' £ ) then v £ | aD = 0. It follows from estimates (I3T251) . (T3T261) . (f3T28l) 
that the set |(Vg, Vg)J is bounded in //'(R) © //'(B) uniformly in e. Therefore the set {(v^, Vg)} 

is weakly compact in H l (R) © // ! (B) and in view of the embedding theorem it is compact in 
L 2 (R) © L 2 (B). Let s' c s be an arbitrary subsequence for which 

v R g //*(R) weakly in H\R) and strongly in L 2 (R), v R | 5D = 

ri^^vu^iii^^. i.,:.r.™ ( 3 - 29 ) 



e=e'— *0 



4 — > vb g // (B) weakly in // (B) and strongly in L 2 (B) 

e=e'->0 



We will prove that 



v = £°f, where v = (v R , v B ) 



(3.30) 



We define the function w £ e //'(D) by the formula 



w £ (x) = (w B (x) - Wr(x)) (f 



\x- x £ i} \-{r-s y - 1 )^ 
^F 1 



+ Wr(x) 



Here wr, wb g C 00 (R' 2 ) are arbitrary functions, supp(w R ) c D, <p : R — > R be a smooth function 
such that (p(p) = 1 as p < 1/2 and <£>(p) = as p > 1. Substituting w £ into (13.271) we get 



I 



(Vv R , Vwr) + v £ R w R - /rWr 



I 



(Vv B , Vw B ) + fey (v B w B - / B w B ) 



rfy + 5(e) = (3.31) 



R 

where 



6(e) 



-I 



(Vv B , Vw B ) + fey (v B w B - /bWb) 



a £ (Vv £ , Vw £ ) + fey(v £ w £ - / £ w £ ) 



It is clear that 



a £ \Vw £ \ 2 dy + ||w £ ||£ 2(GE) < C(e l + e^ 1 ) 



,£||2 
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and due to (13.251) . (13.26b . (13381) we get £ £ -> as £ 0. Taking into account (I3T291) we pass to 
the limit as e = e' — » in (13.311) and obtain 



(Vv R , VvVr) + VrWr - f R W R 



dy + 



(Vv B , Vw b ) + bj (v B w B - /bW B ) 



dy = 



Hence -A£% + Vr = f R and -fcj^va + v B = / B . Therefore (13301) holds. In view of (13301) 
(vr, vb) is independent of the subsequence e' and thus (v R , v B ) converges to (vr, vb) as s —> 0. 
Making the substitution x = + x e . in estimate (13.71) we get 



J a £ |Vv e | 2 rfy + J* |V- 



L 2 (R) 



and therefore in view of ( 13.261 ). ( 13.281) we obtain (recall that y > 3) 

"- e " 2 - 



£->0 



(3.32) 



Taking into account (13391) . (13301) . (13321) we get 

HOT/ - JP£°/fc < K - vrII| 2( r) + K - vbIIL/b^) + 2(l|v £ |li 2(G£) + ||v B |li 2(G£) ) ^ 
and thus C3 is proved. 

Finally let us verify condition C4. Let sup ||/ e ||^ < 00. We denote v E = £, £ f E , it is clear that the 

£ 

set {v £ } £ is bounded in H l (D) uniformly in s. Then the set |(v R , v B )J is bounded in H l (J&) @ H l (B>) 

uniformly in s and therefore the subsequence s' c e and w = (wr, wb) e # ! (R) © ^(B) c < H° 
exist such that 



£=e'-»0 



£=£'-»0 



wr weakly in H l (R) and strongly in H l (R) 
w B weakly in L 2 (B) and strongly in L 2 (B) 



Moreover v E satisfies (13.321) . therefore lim ||X°/ £ _ w \\ 2 = 0. C4 is proved. 

£=£'->0 

We have verified the fulfilment of conditions C\ -C 4 . Thus the eigenvalues /if of the operator £f 
converge to the eigenvalues /i* of the operator X° as e — > 0. But /if ' e (D) = (a^)" -1 - 1, A* = Oui) -1 - 1 
that implies (13.231) . The lemma is proved. □ 



4. Structure of cr(J?l ) 



In this section we prove equality (11.131) . 
At first we show that 

m 

A e \ |J [o-j] <=> ^f(A) e <r0£) (4.1) 

7=1 

— — " d 2 

where <x(J?l) is the spectrum of the operator & = - / a* ? - — — , the function f{X) is defined by 

G3 
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Fig. 2. The graph of the function AT{A) (for m = 3). 



Indeed let A 6 cx(j?l ) \ M cr, . Then there is nonzero F 

7=1 1 J 



that 



fi 

\fmj 



F £ im(J?l - AT, 



e L 2 (R") e_L 2 , P;/(r; (R") such 

7=1 ,m 



(4.2) 



Let us suppose the opposite, i.e. AT(A) i criJK). Then for any g e L 2 (R") there is u e dom(J?l) 
such that 



Jlu - AT(A)u = g 



(4.3) 



ill r 

We set g = f + Y —!—!—. it follows from g3]) that 

^— ' (Tj - A 



;=i 



(7 -AU = F, where t/ 



/ u \ 



(Till + fi 

, uj = (J= \,...,m) 

o-j-A 



We obtain a contradiction with (14.21) . hence AT(A) e cr(j5T). Converse assertion in (14.11) is proved 
similarly. 

It is well-known that cr(Jl) = [0, oo), therefore 



A e o-(Jl ) \ Q {cry} iS Ar(A) > 



(4.4) 



./=! 



At first we study the function AT{X) on R. It is easy to get (see Fig. [2]) that there are the points 
/uj, j = 1 , . . . , m such that 

r(Mj) = 0, j = l,...,m-l 

(Ty < Uj < (Tj+u j = 1, . . . ,m - 1, cr m < yu m < oo 

^ G R \ (J (cry) : AT(A) > 1 = [0, oo) \ \J[o-j,Hj) 
j=i ) V;=l 
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Let us consider the equation AT {A) = a, where a e [0, oo). One the one hand it is equivalent 

/ m r 1 

to the equation I \\icr j - A)\ P m+ \(X) = 0, where P m+ \ is a polynomial of the degree m + 1, and 

therefore in C this equation at most m+ 1 roots. On the other hand on [0, oo) the equation AT{X) = a 

has exactly m + 1 roots (see Fig. [2]). Thus the set {A e C : Af(A) > 0} belong to [0, oo). 

m i i / m \ 

We conclude that A e <t(J1 ) \ U \<Tj\ iff /I e [0, oo) \ I [J [a-j.fiM. Since cr(J?l ) is a closed set 

then the points crj, j = l,m also belong to cr(Ji ). This completes the proof of equality (11.131) . 



5. Proof of Hausdorff convergence 

This section is a main part of the proof: we show that the set cr{J\ E ) converges in the Hausdorff 
sense to the set cr(Jl ) as e — » 0, that is the following conditions dA H [ ) and dB H [ ) hold: 



if /l £ e (r(J?l £ ) and lim /l £ = A then i e cr(J?l ) (A H ) 



e-»0 



for any A e cr(Jl ) there exists A E e cr(Ji E ) such that \imA £ = A (B H ) 

£^0 



5.1. Proof of condition ( A H >• Let A £ e cr(Jl £ ), lim A E = A. We have to prove that A e cr(Jl u ). 

1 ' E->0 



III , . //I , > 

If A e |J |cr 7 j then dA H [ ) holds true since U |cr 7 } c cr(j?l ). Therefore we focus on the case 



7=1 v ' 7=1 

HI 

[<rj\. 

7=1 1 J 

We consider the sequence e w c e, where ejv = ^ N = 1, 2, 3. . . For convenience we preserve 
the same notation e having in mind the sequence s N . 

Taking into account Remark |2~T| we conclude that there exists (f s T" such that A £ e a{^R^ ,£ ). 
We extract a subsequence (still denoted by e) such that 6? £ — » e T". 

£— >() 

Let m £ e dom(J?ly ' £ ) be the eigenfunction corresponding to A £ and such that 

ll" E ||L2,i, E (y) = 1 (and therefore r] £ Y [u £ ] = A £ ) (5.1) 
We introduce the operator IF : H l (Fy) — » H l (Y) such that for each u e H l (F £ ): 



Wu(x) = u(x) for xef 

lin £ w||//i ( y) < ci|m||//i ( ^) 



(5.2) 



It is known (see e.g. EdSl) that such an operator exists. 

Also we introduce the operators II £ : L 2 (U ie j E 5 £ .) —> L 2 (Y) (j = 1, . . . , m) by the formula 

i el £ , xeYf : U £ u(x) = (u) B , 
(recall that Y = \J Y £ ). Using the Cauchy inequality we obtain 

wmw^Y) < q\u\\ l (5.3) 

J 'J 
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It follows from (I5.1I) - (I5.3I) and the embedding theorem that a subsequence (still denoted by e), 
u G H l (Y) and uj e L 2 (Y) (j = 1, . . . , m) exist such that 

TTu E — » u weakly in H l (Y) and strongly in L 2 (Y) 

U E u E —> uj weakly in L 2 (Y) (j = 1, . . . ,m) 

Moreover due to the trace theorem 

U E u E -> u strongly in L 2 (dY) (5.4) 

and therefore u belong to H l e {Y), i.e. 



e->0 

rl, 



Vfc=l,n: w(.<c + e^) = 0ku(x), for x = (xi, x 2 , ■ ■ ■ , 0, . . . , x n ) (5.5) 



T 

k-th place 



JL a2. 



d u 

We denote by Jit the operator which is defined by the operation J{ e Y u = — ) Tt 1 and the 

£d ox k dx t 

definitional domain dom(J?l^) consisting of functions belonging to H 2 (Y) and satisfying ^-periodic 
boundary conditions, i.e. 



Vk = l./i 



w(x + e*) = 6ku(x), 

y + ei) = ft y for * ■ <* ?■ • ■ • ■ *»> 

OX/ OXl t-th place 

/=1 /=1 



It is clear that cr(j?l^) c [0, oo). 
Lemma 5.1. One has 

u e dom(^) and tf^u = AT{A)u (5.6) 
Proof. One has the following integral equality: 

J (a E (x)(Vu E (x), Vw £ (x)) - A E b E (x)u E (x)w E (x)yx = 0, Vw £ e HL(Y) (5.7) 

F 

In order to prove (15.61) we plug into (15.71) a function w e of special type and then pass to the limit as 
e -» 0. 

We introduce some additional notations. Let vj e C 2 (F) (k = 1, . . . , n) be a function that solves 
the problem ( 11.71) in F. We denote by v~ k the function that belongs to C 2 (Y) and coincides with 
Vk in F (such a function exists, see e.g. |fT9l ). Then we extend by periodicity to the whole R" 
preserving the same notation for the extended function. Using a standard regularity theory one can 
easily prove that e C 2 (R"). We set 



k 



v E Ax) = ev k ixE l ) 



Let v E j e C ' E (D E j) (i e Z", j = 1, ... ,m) be the function which is defined in D E . by formula 
(13.101) . supp(v £ ) c D E .. We redefine it by zero in R" \ D E .. Recall that v £ . was constructed in 
Lemma l3~4l as an approximation for the eigenfunction v E . of the operator which corresponds 

to the first eigenvalue Af' E (D E .) and satisfies (13.91) . 
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Let (p : R — » R be a twice-continuously difFerentiable function such that (p(p) = 1 as p < 1/2 
and (p(p) = as p > I. We set 



m i 

tf(*) = l- 



|x - - (r £ - d £ ) 



, x e 



Its clear that 



<p e Ax) = for x e M Bf„ «f (jc) = 1 for x € R" \ M (Gf ; U Bf ,.) 



7=1 



\y=i 



(5.8) 



|D>f| < Ce"^ (M = 0, 1, 2, . . . ) in \J Gf 

7=1 



We cover R" by the cubes 

Yf = [x 6 R" : < x A < (i + l)e + e 3/2 } , t e Z" 

Let (x)} . be a partition of unity associated with this covering, that is 

<Af(x) e C 2 (R"), <#(*)<!, = 1, 0f(jc)= lifxe yf\[Jyf, ^(x) = 0ifxgFf 

r'eZ" fri 

Moreover, analyzing a standard procedure of the construction of the partition of unity (see e.g. Il2~0l0 
we can easily construct the partition of unity satisfying the following additional conditions 



Vi 6 IT, Vx e R" : tf(x) = i/f Q (x + is) 



£>> £ (x) < Cs~ 3a/2 (a = 0, 1, 2) for x e Y n 
We consider the function w E of the following form: 



U 7 < 



\x) = J]tf(x)(g E (x) + h E (x)) 



(5.9) 
(5.10) 

(5.11) 



where 



g%x) = g(x" E ) + tf(x) I J] p-(xf' s )(x k - 4 s - v £ (x))+ 



.k=l 



dxu 



+ 



1 " 



d 2 g 

2 ^ dxfcd*/ 



- xf - v £ (x))(x, - x Ue - v E (x)) 



fcf(x) = X (hjtfj) - gix^j 

7=1 

Here x ! ' £ is the center of Y E , g(x), hj(x) are arbitrary functions from C 2 (R") satisfying 
VxeR", V/ = (i u ...,in) gZ" : 



(5.12) 
(5.13) 



+ = 0? ■...•#£(*) 



+ i) = 0'j 1 



0£fy(x), ; = l,...,m 



(5.14) 



Remark that = on dG E t . Taking this into account we conclude that w £ (x) belongs to C 2e (R") 
on 1 



and in view of (15.9k (15.141) and the periodicity of v\ we get 

Vx 6 R", Vi e Z" : w £ (x + i) = flf • . . . • 0j> £ (x) 
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We also introduce the notations 

g £ = J] ^(x)gf(x), h £ = £ 

(eZ" jeZ" 

The function w £ belong to Hi(Y). In order to obtain the function from HL(Y) we modify w £ 
multiplying it by the function which is very close to 1 in Y as s — > 0. Namely, we define the 
function 1 £ e C°°(K") by the following recurrent formulae: 

1 (x\, . . . , x n ) = A n {x\, . . . , x n -i)x n + B n (x\, . . . , x n -i), 



a = 2, . . . ,n 



It is easy to see that 



Finally we set 



I B a (X\, . . . , X a -\) - Aq._i(Xi, . . . , X a -2) x a-l + B a -i(Xi, . . . , X a -2), 

\A a (xi, . . . , x a -i) = {6 E a l6 a - \)B a {x\, . . . , x a -i), 
B x = 1, A\ = WjWi - 1. 



max |1 £ (jc) - 1| + max |Vl £ (x)| -> 
1 £ e #1/10, where # £ /# := (# £ M, . . . , 6 E je n ) 



(5.15) 



w £ (x) = w £ (jc) + (l £ (x) - l)w £ (x) 

It is clear that w £ e HL(Y). 

8 s 

Substituting w £ into (15.71) and integrating by parts we obtain 



- A £ u £ w £ )b £ dx+ 



+ 



J u E ^-ds + J (a £ {Vu £ , V((l £ - l)w £ )) - A £ b £ u £ {\ £ - l)w £ )jx = (5.16) 



BY 



Further we will prove that the second and the third integrals in (15.161) tend to zero as s — > 0. 
Now we focus on the first integral in (15.161) . Using Lemma U31 and taking into account (15.11 ), (15.81 ) 
we obtain the estimates 

I WIlLtpj) * ^ (5.18) 
Since = in G £ . then in view of (fO) . (IBTTTl) . (T5TT81) 

|( M s ,^ s w%^ (G£y) |<C e (5.19) 

Similarly we obtain 

lim( M £ ,w £ ) L2iE(G , ) = (5.20) 

e— >0 ' 7 

We denote 

Ff = {jc G F £ : is + e 3/2 < jc t < (i + l)e) , F £ = [J F £ 



It is clear that Ff = Ff \ I jj Yf ). 
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Firstly we estimate g E in F £ \ F £ . We represent g E in F E \ F E Y in the form 



k=\ \ 1=1 



g%x) = J] U*) ft* + 2 - xf) + \ E - *£•)(*, - 4 £ ) - g(x) 

~ E I ^ + E " ^' e) - I + ? - s,««K(*K (*) + 

/ JU=1 

+ - X + 5 E g,«(*K(*K(*) (5.2i) 



k=\ 



k,l=\ 



ds d 2 s 
Here g'- E = g(x'' E ), g, k (x) = -£-(*), g'f = g*(x J ' £ ), = T-f-W, ^ = g e H (^ e ). It follows 

OXk ' ' OXkOXi 



from (Ell) , (E2p that |Ag £ (x)| <Cforxe^\ F £ . Since dist ( Q £>?., d7f ) > ks then /i £ = in 



Fy \ Fy when e is small enough and therefore 



F*| < C 2 s 



1/4 



Similarly we obtain 



\im(u\w s )^ (FeA FS) 



Let us study g E and /z £ in F E . It is clear that 

n 

A/ = ^ g^(V(x* - v £ ), V(x, - vf)) for x € Ff 



(5.22) 
(5.23) 

(5.24) 



In view of Lemma I3TT1 and the Poincare inequality one has the following estimate: 



u E - (u E )j 



_ +s" 

L 2 (Ff) 



{Il E u E ) Y c - (u E h 



+ \\U E u E - (U E u E ) Yf \\l 2(YD < Ce 2 \\Vn E u E \\ 2 L2{Yf) (5.25) 



Using (15.241) . (15.251) and the Poincare inequality we get 



-I 



j (V(**-v k ), 



V(x/ - v,))dx 



V F 



=-e^ f i r 

kj=l X 



d 2 g 



■dx + o{\) 



iel £ 

kj=i i 



+ o(l) = 
d 2 g 



dx k dxi e^O Z-J J dx^dx/ 



Jjc (5.26) 



In the same way using Lemma 1X21 (for v E = U E u E ) we obtain 

^^S £ ) L2 ^ y) = YjS^ £ ){u £ )FAF\s n + o(l) -> \F\ f M (x)g(x)Jx (5.27) 

(here we also use the estimate £"|Kw £ )™ - (w £ )HI 2 < C£ 2 ||Vn £ w £ || 2 lve . which follows from Lemma 



ED- 
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Let us study h E in F y . Integrating by parts and taking into account the form of the functio n vf . 
(in particular, we have the estimate ||J?T e v^.||L 2 (yp < Ce n ), the Poincare inequality and Lemma [3T21 
we obtain 

m 

+ Z (" £ - <«v ^L,^) = Z ^ M Z< m > - w) «" + °(d £ 

/eJ £ ' ' j=l ieI E 

m p 

^aj\dBj\ I u(x)(g(x)-hj(x))dx (5.28) 



e-»0 

7=1 



(here r = |x - jc?.|). In the same way we get 

Urn( M s ^% 6e( ^ = (5.29) 

Let us study h E in B Y (g E = in 5^). Integrating by parts and using the Poincare inequality we 
obtain 

iu E , wm^m) = - Z Z ( " e H -j 1 (W) - Js + °( 1 ) = 

7=1 'eJ £ ^ r 

^aj\dBj\ I n^ e (x)(^(x) -,?(x))Jx + o(l) -» ^aj\dBj\ I h ; (x) (/j/jc) - g(x)) dx (5.30) 



7=1 y 7=1 

In the same way we get 



m p 

X\m(u\h £ ) L2 ^ = J^lBjlbj I Uj(x)hj(x)dx (5.31) 
J= l a 

Finally, let us estimate the remaining integrals in (15.161) . One can easily obtain that 

T] £ Y [w E ] + \\w% bE(Y) < C 

and therefore in view of (15.151) 

lim J (a E (Vu E , V((l £ - l)w e )) - A E b E u E {l E - l)w e )jx = (5.32) 

Y 

dw E 

It is easy to see that the function p E = is bounded in L 2 (dY) uniformly in e and therefore 

on dY 

there is a subsequence (still denoted by s) and p e L^idY) such that 

p E -> p weakly in L 2 (dY) (5.33) 

Moreover it is clear that Vk = l,n: p E ix + et) = -(F k p E {x) for x = (x\, X2, . . . , 0, . . . , x n ). Therefore 



T 

k-th place 



Vk=l,n: p{x + et) = -6kp(x) for x = (x\, xi, . . . , 0, . . . , x„) (5.34) 



T 

A--th place 
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Taking into account (15741) . (1531) . (15.331) . (15.341) we get 



r dw £ r 

I w £ — — ds = \ up ds = (5.35) 
J on J 



lim 

£->o J on 

dY dY 



Then taking into account (I5TT9I) . (157201) . (157221) . (157231) . (15726l) - (15732l) . (157351) we pass to the limit 
in (15.161) and obtain the equality 

-u{x)\F\ Y a*'— f-(x) - A\F\u{x)g{x) + Y L|d5 7 -|(g(x) - ^(x))«(x) + 
V W=l ax k° x l j=l v 

+a 7 |35 ; |(/i ; (x) - g(x))uj(x) - A\Bj\bjUj(x)hj(x)Jj dx = (5.36) 

Recall that g, hj e C 2 (R") are arbitrary functions satisfying (15.141) . 

Plugging g = 0, /i ; = for j # into (15.361) and taking into account the equality \dBj\ = \Bj\nr~ l 
we get 

o~k 

lit = u, k &{!,..., m} (5.37) 

cr k -A 

Then setting hj = for all j = 1 , . . . , m, integrating by parts and taking into account (15.371) we get 



( n „ „ A 

dx = (5.38) 



where the function T(A) is defined by (11.91) . 

Equality (15.381) is valid for an arbitrary g belonging to C°°(R") and satisfying (15.141) . It is clear 
that the set of such functions is dense in H-(Y). Therefore equality ( 15.381) implies (15 .6b . Lemma 

■ 9 

D.ll is proved. □ 
Lemma 5.2. u ± 0. 

Proof. Let us introduce the spherical coordinates (r, 0) in D? . and the function by the formula 

u\fp, 0) = <M £ ) sf;(P ), where Sf/p) = {x e R" : |x - x £ .| = p) 
One has the following Poincare inequality: 

iaS VqK = ' 

by p from to r £ - rf £ and summing by i we get 

Yj ^ - M PW> * ^ 2 ||V M £ ||i 2(UBf) < C l£ 2 i £ (5.39) 



- C P 2 H V © M£ |lL2(Sf ; (p)) - C l e2 |l Vw£ |lL 2 (Sf J (p)) 

(here V© is a gradient on S £ (p): for example in the case n = 2 one has V@m = Trf|). Integrating it 



We denote u £ . = wf. - {u E ) s ^. Clearly uf. e dom(J?fc £ ) and 
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Recall that A g Q Then in view of Lemmas I3.4L [331 /l £ £ cr(j?T ' £ ) when e is small enough. 
Therefore we have the following expansion: 



uf ; = J] /5(e), where 7* (e) = v? (£>*■)- 



(i £ <w £ ) 5 ,,vf(D £ .)) 



(^ £ (D £ ) - A") 



(5.40) 



Here {vf (£> £ )] is a system of eigenfunctions of corresponding to {/if £ (.D £ .)) _ and such 
that (v?(PfX vf(Df .)) =0ifk±l. 

Using Lemmas I3T21 13.51 we get (for j e {1, . . . ,m}) 



CO 



e" -> 



(5.41) 



As in Lemma I3~4l we denote v £ . = vf(D £ .) assuming that v £ . is normalized by condition ( 13.91) . 



Using estimates (13.161) . (13.181) and Lemma I3T21 we get 



< M £ >< 



(10 



(5.42) 



as e -» 0. It follows from d5.40h - d5.42h that 
Similarly we obtain 



^ 2 |5;II|M" 2 



(5.43) 



r (u £ ) sl A\Bj\ 
I u? Jjc : — : — e n as s -> 



try -/I 



Using (|3A6b . (I318T) . (15431) . (f5~44|) and Lemma[32lwe get 



(5.44) 



Zll« 



e->0 



A 2 \Bj\ 2A\Bj 



"vibe*.) + 2 <" £ >s 



(x)dx + 



+ 



(o-j-A) 2 o~ j — A 



+ \Bj\ 



(u £ h 



cr, 



■ \B £ \s" 



s->0 



o-j-A 



\L 2 (Y) 



Using the Poincare inequality and Lemma I3T21 one can easily prove that 



^%, b ^ y) = \F\J](u E h f e n + o(l)^ Q \F\. 



\\L 2 (Y) 



Furthermore in view of Lemma 1331 



lim||w £ || 2 rrE , = 



(5.45) 



(5.46) 



(5.47) 
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Finally taking into account ( 15.391 ), d5.45h - d5.47b we obtain 



1 = \\u 



e\\2 



\L 2M s(Y) 



mi \2 



and therefore u ± 0. Lemma [521 is proved. 



□ 



It follows from Lemmas |5 . 1 H5 . 21 that AT(A) belong to the spectrum cr(M e Y ) of the operator tR e Y . 
Therefore AT (A) e [0, oo) and in view of (1441) A e cr(Ji ). Condition dA^j ) is proved. 



5.2. Proof of condition ( |B7J ). Let /I g <x(j?l ). Let us prove that there is A E e cr(^l e ) such that 
\imA £ = A. 

£->0 

We assume the opposite: the subsequence (still denoted by s) and 6 > exist such that 



dist(i, cr(J?l £ )) > 6. 



(5.48) 



Since /I e cr(Jl ) then the function i 7 



\fmJ 



G L 2 (R") e_ L 2 , Pj/crj (R") exists such that 

j=l,m 



F i - AT), where I is the identical operator (5.49) 

It follows from (15481) that A e R \ cr(j?l e ). Then im(J?l £ - ,11) = L2 >fce (R") and hence for an 
arbitrary / e e L 2ifo£ (R") there is the unique u £ e dom(J?l e ) such that 

Jl £ u £ - Au £ = f (5.50) 

We substitute the following f £ {x) e L 2j ^(R") into (15.501) : 



fix) 



{f) Yt , x e Ff, 
(fj) Y c, xeB £ r 



0, 



It is clear that the norms H/^Il^cr") are bounded uniformly in e. Then in view of (15.48b u £ satisfies 
the inequality 

Furthermore 



\\Vu 



£||2 

L 2 (R") 



<ll/ £ 



l^llz^CR-) + \A\ • \\U £ \\1 



2,6* (R") 



< C 



Hence a subsequence (still denoted by s) and u e H l (R"), Uj e L^R") such that 

n £ « e — » u weakly in H l (R") and strongly in L 2 (G) for any compact set G c K" 
U £ u £ uj weakly in L 2 (R") (j = 1, . . . , m) 

where II s , ITj (j = 1, . . . , m) are the operators introduced above in the proof of condition dA H b - 

o 

For an arbitrary w £ e C M (R") one has the following integral equality: 

J (a WW , V») - - OWWU))^ = (5.51) 
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We substitute into (15.511) the function w E of the form (15.1 ll) - (15.13i but with g, hj e C°°(R"). Making 
the same calculations as in the proof of condition ( |A H | ) we obtain 

- u{x) £ ahIFI^^C*) - A\F\u(x)g(x) - \F\f(x)g(x)+ 



k,i=i 



+ ]T l aj \dBj\(g(x) - hj(x))u(x)+ 
j\dBj\(hjix)-g(x))uj(x) - A\Bj\bjUj(x)hj(x) - \Bj\bjfj(x)hj(x)\ 



dx = (5.52) 



for an arbitrary g, hj e C^IR' 1 ). It follows from (15321) that 



U 



"i 



6 dom(J?l ) and fPu - AU = F 



We obtain a contradiction with (I5.49I ). Condition ( |Bh| ) is proved. 

6. End of proof of Theorem [TTI 

In general the Hausdorff convergence of cr(Jl E ) to cr(Jl ) does not imply ([TTTO - ([TTTTlH . However 
if we prove that cr(Jl E ) has at most m gaps in [0, L] when s is less some s L then this implication 
holds true. More precisely the following simple proposition is valid. 



Proposition 6.1. LetS E = [0,L] \ U(o^,^) , S = [0,L] \ U L < oo and 



\j=i 

< off, a E < /3 E < a E j+l , j= l,m s -l, < e < L 

< ori, ar ; - < /Jy < ary+i, j = l,m - 1, a m < L 
m £ < m 

S e converges to S in the Hausdorff sense as e — » 
77zen m E = m when s is small enough and 

Vj = 1, . . . , m : lima^ = a,-, lim^ = /3,- 



7=1 



We introduce the notation [a7(s), aUe)] := U f4' e ( y o)|- 

6eT" 



Lemma 6.1. \ima + ^(e) = oo 



e->0 



iVoo/ In the same way as in the proof Lemma [331 we obtain the following equality 

\ime 2 A N k E {Y E ) = A k , k= 1,2,3... 



(6.1) 



For example, the set a* :- cr(Jl ) n ( (J e(£ + ^)1 I also converges to cr(j?l ) in the Hausdorff sense, but the 



number of gaps in cH 1 n [0, L] tends to infinity as e — > 0. 
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where [A k } km are the eigenvalues of the operator A which acts in the space L 2 (F) ®_L 2>b .(Bj) and 

j=i,m 

is defined by the operation 



F 



A = 







b~ l A N 












(here and A^ are the Neumann Laplacians in F and Bj). It is clear that Aj = for j = 
l,...,m+ 1 while A m+2 > 0. Then using (16.11) and taking into account (12.31) we get 



^^ m+2 (s) > lirrU +2 <T e ) = A m+2 lims- 



OO 



Suppose that there is a subsequence (still denoted by s) such that the numbers a* n+l {e) are 
bounded uniformly in e. Let the numbers L, L\ be such that fj, m < L <L\ and < L. Since 

lima" 2 (e) = oo then <3~ Js) > L\ when e is small enough. Hence cr(Jl E ) n [L, = when e is 



small enough. We obtain a contradiction with condition ( |B H | ) of the Hausdorff convergence. Thus 
lir n< +1 (£) = oo. □ 

e— *0 



Lemma [6J] implies that for an arbitrary L > the spectrum <x(j?l £ ) has at most m gaps in the 
interval [0, L] when e is small enough: 



cr(Jl E )n[0,L] = [0,L]\[j( ( T E j ,M E j) 

where (cr £ ,//p c [0, L] are some pairwise disjoint intervals, rrf < m. Here the intervals are 
renumbered in the increasing order. 



We have proved that <j(J{ E ) converges to [0, oo) \ (J (cry,//,) in the Hausdorff sense as s — » 0. 



./=! 



Let L be an arbitrary number such that L > /j, m . Then, evidently, cr(Jl E ) D [0, L] converges to 



[0, L] \ ^ U (o-j,jUj) ) in the Hausdorff sense. By Proposition ^. ll m £ = m when e is small enough and 



Vj = 1, . . .,m : 



lim cr E 

s->0 7 



0> lim^ = 



Theorem ll.ll is proved. 



7. Proof of Theorem 1 1.21 
Substituting a h bj (11.121) into (fL~8T) we get 



(i.e. the first equality in (10.51) holds) and 



i=l,m|i# j 



(7.1) 
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Recall that jij (j = l,m) are the roots of equation (11.91) , therefore in order to prove the equalities 
/iy = fij ( j = 1 , m) we have to show that 

m 

n=l,...,m: S\-£l— = \ (7.2) 



Let us consider (17.21) as a system of m linear algebraic equations (pj, j = 1, . . . , m are unknowns). 



It is clear that (17.21) follows from the following 

Lemma 7.1. 77ze system ( 17.21 ) has the unique solution p u . . . ,p m which is defined by ( 17. 

Proof. We prove the lemma by induction. For m = 1 its validity is obvious. Suppose that we have 
proved it for m = N - 1 . Let us prove it for m = N. 

Multiplying the k-th equation in (17.21) (k = 1, . . . , N) by - and then subtracting the N-ih. 
equation from the first iV — 1 equations we obtain a new system 

N-l 

V* = l,...,tf-1: V-^— = 1 

where the new variables pj, j = 1 , . . . , N - 1 are expressed in terms of pj by the formula 

a m — a ; 

Pj:=Pj- -, ; = l,...,iV-l (7.3) 

J3 N ~ ocj 



Hence pj, j = l,N - 1 satisfy the system (17.21 ) with m = N — 1. By the induction 



(•=1,AT-1|(^ 



It follows from (O, G4J) that p ; - (j = 1, . . . ,N - 1) satisfy (O) (with m = N). The validity of 
this formula for p N follows easily from the symmetry of system (17.21) . Lemma FTTH is proved. □ 

Theorem [L2]is proved. 
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